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Abstract 

We compute various form factors involved in i? — D^*^ transitions based 
on the perturbative QCD formalism, which includes Sudakov effects from the 
resummation of large radiative corrections in a heavy-light system. A two- 
parameter model wave function for D^*^ mesons is fixed using data of the 
nonleptonic decays B D^*^7r, from which the ratio of the decay constants 
Id*/ fn = 0.92 is obtained. We then derive the spectrum of the semileptonic 
decay B D*ii> in the fast recoil region of the D* meson, and extract 
the CKM matrix element \Vcb\ = 0.043 x (0.12 GeV/Zfi) x (0.14 GeV/fn), 
Jb and fo being the B and D meson decay constants, respectively. Here 
we adopt the convention with the pion decay constant = 93 MeV. With 
these outcomes, we evaluate the decay rate of B ^ DDg, and estimate the 
ratio fos/fD = 0.98 from data. Contributions of internal W^-emission and 
IV-exchange diagrams are briefly discussed. 



1 Introduction 



Recently, the perturbative QCD (PQCD) formalism including Sudakov ef- 
fects has been shown to be applicable to heavy meson decays, which were 
usually regarded as being dominated by nonperturbative dynamics. The 
breakthrough is attributed to the resummation of large radiative corrections 
in heavy-light systems such as a -B meson containing a light valence quark. 
This resummation, which was first performed in [|l| for the semileptonic de- 
cays B — > 7r(p)£z/, improves the applicability of PQCD to these heavy-to-light 
transitions. It was found that PQCD predictions are reliable for the energy 
fraction of the pion above 0.3, from which the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element \Vub\ can be extracted, once experimental data are 
available. A similar formalism was applied to B ^ D decays in the fast recoil 
region of the D meson p| , and found to be self-consistent for velocity transfer 
above 1.3. The resummation technique was further employed to study the 
inclusive decays B and B X„£z/ 0, and the Sudakov effects at the 

end points of spectra were examined @]. 

In this paper we shall extend the analysis in 0, and evaluate all the 
form factors involved in 5 ^ D'^*'> transitions at the high end of velocity 
transfer rj. The B meson wave function is determined by the relativistic 
constituent quark model For fast D^*^ mesons, a convincing model wave 
function is still not yet obtained. We propose a two-parameter D'^*'> meson 
wave function, and fix the parameters using data in the large rj region from 
the nonleptonic decays B D^*^7i One of the parameters corresponds 
to the normalization constant, and the other controls the shape. With these 
phenomenological inputs, our analysis is free of the ambiguity from nonper- 
turbative effects. We compare our results to those from other theoretical 
approaches, for example, heavy quark symmetry (HQS) combined with 
0{as) and 0{1/M) corrections in ||^ and overlap integrals of heavy meson 
wave functions in [10|. 

We derive the spectrum of the semileptonic decay B D*lv from the 
transition form factors. The CKM matrix element \Vcb\ is then obtained by 
fitting our predictions to data [O at large where PQCD is reliable. Our 
approach to the extraction of |V"cb| differs from those in the literature, where 
the behaviors of the form factors at 77 = 1 are employed. The relevant nonva- 
nishing form factors take the same functional form ^(r^), which is normalized 
to unity at zero recoil, ^(r^ = 1) = 1, because of HQS. The quantity |V^f,|^(l) 
can be extracted from experimental data, if the behavior of ^ above zero re- 
coil is known. However, ^ is thought of as being uncalculable in perturbation 
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theory, and thus the extraction depends on how to extrapolate ^ from r) = 1 
to 7] > 1. Hence, different models of ^(77) lead to different values of \Vcb\ 
[§, |TI|. We argue that ^ is in fact calculable in the fast recoil region. The 
nonperturbative wave functions, fixed by data of other decay modes, provide 
model-independent extraction of \Vcb\- 

The consistency of \Vcb\ determined in this paper with currently accepted 
values justifies the PQCD analysis of B meson decays, especially B ^ ti 
decays. Predictions for the spectra of the decays B — > 7r(p)£i/ in are 
then convincing, and can be used to extract \Vub\- On the other hand, HQS 
requires only the normalization of heavy-to-heavy transition form factors at 
zero recoil. The PQCD formalism, however, gives information near the high 
end of 7]. Therefore, these two approaches complement each other. Fur- 
thermore, our formalism is applicable to the evaluation of nonfactorizable 
contributions from VT-exchange diagrams, which remains a challenging sub- 
ject in the study of nonleptonic B meson decays. 

We develop the factorization formulas for all the B D^*^ transition 
form factors in Sect. 2. In Sect. 3 the two-parameter D^*^ meson wave 
function is determined by fitting the data of the decays B D^*^7r. The 
B D*Cv spectrum is derived, and \ Vcb\ is extracted from experimental data 
We show in Sect. 4 that nonfactorizable contributions from VT-exchange 
diagrams are negligible. Sect. 5 is the conclusion. 



2 Factorization Formulas 

The six form factors ^j, i = +, — , V, Ai, A2 and ^3, involved in i? D^*^ 
transitions are defined by the following matrix elements: 

{D{P2)\v^\B{p,)) = ^MBMoii+mvi + + i^mvi - V2r) , 

{D\P2)\A^\B{Pr)) = ^MsMoAUMir] + - UivV ■ 

-eA3(r/)e*-t;i<]. (1) 

The momentum Pi {P2), the mass Mb {Mjj(y)) and the velocity Vi (^2) of 
the B (D^*^) meson are related by Pi = MbVi (P2 = ^d(*)'^2)- The velocity 
transfer 77 = f 1 ■ f 2 has been introduced before, whose expression in terms of 
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the momentum transfer = {Pi — P2Y is given by 



^ Ml + Ml,,, - 

e* is the polarization vector of the D* meson, satisfying e* ■ V2 = 0. The 

vector current and the axial vector current are defined by = cj^b 
and A^^ = c'-f'^'j^b, respectively. 

In the infinite mass limit of Mb and Mj^(y) , the six form factors have the 
relations 

e+ = = u = Us = e, u = u = 0, (3) 

where ^ is the Isgur-Wise (IW) function |^ mentioned in the Introduction. 1^ 
is normalized to unity at zero recoil from HQS. For the behavior of ^ above 
zero recoil, there is only model estimation from the overlap integrals of heavy 
meson wave functions [§, p!0[] . 

We work in the rest frame of the B meson, in which Pi is written, using 
light-cone components, as Pi = Mb/ V^{1, 1, 0^), and P2 has the nonvanish- 
ing components 



= '-^^M,,., . (4) 

AsT] 1 with P^ = P2 = Mb,(,) /a/2, the D^*'^ meson behaves like a heavy 
meson. However, in the large rj limit with P2 S> M^{,) / a/2 ^ P2 ■, the D^*^ 
meson can be regarded as being light 0. We argue that ^ is dominated 
by soft contributions in the slow D^*) meson limit, where the heavy meson 
wave functions strongly overlap, and factorization theorems fail. However, 
when the D^*-* meson recoils fast, carrying energy much greater than Mbi(,), 
B D^*^ transitions are then similar to -B — >■ vr ones ||^ as stated above, 
and PQCD is expected to be applicable. In this paper we shall show that 
the PQCD formalism including Sudakov effects gives reliable predictions for 
^ in the large 77 region. 

The PQCD factorization formulas for the B ^ D transition form factors 
have been derived in 0. Here we summarize the idea. The factorization of 
the matrix elements in eq. (0) into the convolution of a hard scattering am- 
plitude with the B and D*-*^ meson wave functions is shown in fig. la, where 
the b and c quarks are represented by the thicker and thick lines, respectively. 
ki (^2) is the momentum of the light valence quark in the B (D^*^) meson. 
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satisfying A;^ ^ ~ 0). ki has a large minus component fcf , which defines 
the momentum fraction Xi = /P{~ , and small transverse components ki^, 
which serve as the infrared cutoff of loop integrals for radiative corrections. 
Similarly, k2 has a large component A;^, defining X2 = k2 /P^, and small \i.2T- 
A simple investigation shows that large logarithms arise from radiative 
corrections to the above factorization picture. In particular, double (leading) 
logarithms occur in the reducible corrections illustrated by the 0{as) dia- 
grams in fig. lb [m , when they are evaluated in axial gauge. In order to have 
a reliable PQCD analysis, these double logarithms must be organized using 
the resummation technique, which leads to the evolution of the B (D^*^) 
meson wave function 0^ {<Pd(*)) {^t)- We quote the results as follows 



')b{xi,PiM,Ij) ~ 0B(a;i)exp 



-s{k-,M)-'2 r %{a,m 

Jl/bi u 



0D(*)(^2,^2,&2,/i) ~ 0£,(*)(x2)exp 



■/^ C?/i 

/bi fi 

-s{k+,h)-siP+-kt,h) 



-2 r —iMfi)) 

Jl/b2 fi 



(5) 



where bi (62) is the Fourier conjugate variable to kiT {k2T), and can be 
regarded as the spatial extent of the B (D^*^) meson, /x is the renormalization 
and factorization scale. The exponent s organizes the double logarithms 
from the overlap of collinear and soft divergences, and the integral, with 
the quark anomalous dimension 7 = — a^/vr as the integrand, groups the 
remaining single ultraviolet logarithms in fig. lb. The expression of s is very 
complicated, and exhibited in the Appendix. 

In the resummation procedures the B meson is treated as a heavy-light 
system. The D^*^ meson is, however, treated as a light-light system as in- 
dicated in eq. (H), because we concentrate on the fast recoil region The 
initial conditions of the evolution, i = B, D and D*, are of nonpertur- 
bative origin, satisfying the normalization 

1 f. 

(f)i{x)dx = — , (6) 



2^3 

with fi the meson decay constants. The initial condition in eq. should 
be written as 6, 1/6), and (j){x) is in fact an approximation. We have 
neglected the intrinsic transverse momentum dependence denoted by the 
argument b, and PQCD corrections proportional to as{l/b), because these 
two effects cancel partially M. 
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The evolution of the hard scattering amphtude H from the summation 
of single ultraviolet logarithms is expressed as 



H{k^, kt, h, 62, /i) ~ H^'^HK, bi, 62, t) exp 



-4 



(7) 



where the variable t denotes the largest mass scale of H. We have approx- 
imated H by the 0{as) expression H^^\ which makes sense if perturbative 
contributions indeed dominate. Combining the evolution in the above for- 
mula, we obtain the complete Sudakov factor e""^, where the exponent S is 
given by 



S{k^ ,/;;+, 61, 62) 



s{k:[,bi) +s{k+,b2) +s{P^ 
1 



kt,h 



2 ' "2 J 



Pi 



-ln(6iA) -ln(62A) 



(8) 



with (3i = (33 — 2nj)/12, = 4 being the flavor number. The QCD scale 
A = Aqcd will be set to 0.2 GeV below. It is easy to show that e""^ falls 
off quickly in the large b, or long-distance, region, giving so-called Sudakov 
suppression. 

With the above brief discussion, the only thing left is to compute H^^^ for 
each form factor. The calculation of H^^'^ for ^± has been performed in 
In a similar way we derive H^^^ for other ^'s. The factorization formulas for 
all the transition form factors in b space, with Sudakov suppression included, 
are listed below: 



^4 



IGnCpy MbMd J dxidx2 J bidbib2db2(f)B{xi)(f)D{x2) 

x[{Mb + X2C+MD)h{xi,X2MM) + {Md + XiC+MB)h{x2,XiMM)\ 

X exp[-S'(A;f , A;^,6i,62)] 
-IQtiCfJ MbMd ' 



(9) 

dxidx2 / bidbib2db2(f)Bixi)(l)£){x2) 
Jo 

xC_[x2Md/i(xi,X2,&1,&2) " XiMb/i(x2, Xi, 62, &l)] 

xexp[-Sik^,kt,bub2)], (10) 
1QtiCf\J MbMd* j dxidx2 J &i(i6i62'^&2 0_B(a;i)0D* (2^2) 

X [{Mb - X2CiMD')h{xi, X2, bi, 62) + (M^* + xiC2Mb)/i(x2, Xi, 62, &i)] 
xexp[-Sik^,k+,bi,b2)] , (11) 

16tcCf\/ MbMd* J dxidx2 J 6i(i6i62'^&2 0b(xi)0d* (^2) 
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x[(Mb - X2(3MD*)h{xi,X2,bi,b2) + {Md* + Xi(iMB)h{x2,Xi,b2,bi)] 

(12) 



X exp[-S{k^ ,k^,bi,b2)] , 

= -IQttCfJ MbMd* [ dxidx2 [ bidbib2db2<f)B{xi)4>D*{x2) 
* Jo Jo 

xxi(5MBh{x2, xi, 62, bi) exp[-S{k^, k^, bi, 62)] , 
with the constants 

' 1 
2 

3 
2 



(13) 
(14) 



C- 
Ci 

C2 

Ca 



3 

^"2 + 



1 
2 

r]-2 



V h 

' 2 



Iri - 1 
r7 + 1 
+ 1 

7] — 1 



2 2v/^^^^ 
1 



2vV^ 

2 — rj — y/f]'^ — 1 

ri + l 



2{r] + l) 



(15) 



Ci? = 4/3 is the color factor. The function h, coming from the Fourier 
transform of H'^^\ is given by 



h{xi,X2MM) = as{t)Ko {JxiX2CMBMj:)Mbi 



X 



e{b^ - b2)Ko Jx2CMbMi,(,M h Vx2CMbM^(.)6: 



+e{b2 - bi)Ko i^^X2CMBMj,Mb2j lo [^JX2CMBMJ,^.M 

(16) 

with the constant C = + ^ 1- The scale t is chosen as [0 

(17) 



t = ma.x{^XiX2CMBMD{*)A/bi, I/62) • 

Note the equality of and ^^3- We argue that C,v will differ from ^^3, if 
higher-order corrections to the initial condition H{ki, k2, &i, &25 1) are taken 
into account. 
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3 Determination of \Vcb 



Before proceeding to the numerical analysis of eqs. (^)-(|T^), we discuss how 
to fix the B and D^*^ meson wave functions. For the B meson, we adopt 
the wave function from the relativistic constituent quark model 0, (pB^x) = 
/ci2kT/(47r)VB(a;,kr), with 



^Bix,kT)=NB 



Cb + 



Ml 



+ 



1 -2 



X 



x] 



(18) 



The normalization constant Nb and the shape parameter Cb are determined 
by the conditions 



dx 



dx 



d'^k 



(47r)2 

T r i 



(An) 



')B 



^(a;, ky) 
{xM)f 



f. 



B 



2^3 
1 

2 ■ 



(19) 



We obtain Nb = 49.5 GeV^ and Cb = -27.699845 GeV^ for Mb = 5.28 
GeV [|12| and the B meson decay constant Jb = 0.12 GeV from the lattice 
calculation [r^|. The B meson wave function is then given by 



^b(x 



B 



x(l 



X) 



167iMI + Cb{1 



x) 



(20) 



For the D^*^ meson, we propose the following model, which possesses the 
same functional form as eq. (pOD, 



AT^t.) x(l - xy 



x] 



(21) 



Equation (|T9|) is not appropriate for the determination of Nj^m and Cj^i*), 
when the D^*'' meson recoils fast 0. Hence, we shall take an alternative ap- 
proach. It is known from PQCD factorization theorems that wave functions 
are universal, or process-independent. It hints that we can fix the parame- 
ters using data of any B D'^*^ decays, such as the two-body nonleptonic 
decays B D^*^7i, for which factorization theorems should work best. With 
these phenomenological inputs at specific kinematic points rj = r^max, we then 
predict the behaviors of all ^'s in a finite range of t]. 

We assume the vertex factorization hypothesis for the following analysis, 
which has been shown to be consistent with current experimental data [|1^ . 
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We list the formulas for the decay rates of various B D^*^ decay modes 
involved in our study. They are 



327r 



2r 



1 + r 
1 — r 



327r 



1 + r 
1 — r 



(2r)3 

^Ai(^?max) - (r^A2(^?max) +^A3(^max)) 



(22) 



for the nonleptonic decays, with the constants r = Mjj(*)/Mb and 7].^ 
(l + r2)/(2r), and 



X <^ 2(1 - 27]r + 



(23) 



for the spectrum of the semileptonic decay D*^i^u p|. For nonlep- 

tonic decays, there exist additional important corrections from final-state 
interactions with soft gluons attaching the D^*^ meson and the pion. It has 
been argued that these corrections produce only single logarithms, which 
cancel asymptotically 0, and are thus not considered here. 

We adopt Gp = 1.16639 x 10~^ GeV~^ for the Fermi coupling constant, 
\Vud\ = 0.974 for the CKM matrix element, Md = 1.87 GeV and Md* = 2.01 
GeV for the D and D* meson masses [0, respectively, r^o = 1.53 ps for 
the B^ meson lifetime [rT|, f^^ = 93 MeV for the pion decay constant and 
fi) = 0.14 GeV for the D meson decay constant from the lattice calculation 
jl^ . The four parameters we shall determine are Cu {Nd is fixed by fn from 
eq. (D), Cf)*, fjj* (or Nd* equivalently) and At the same time, we 

have four constraints from experimental data: the branching ratios B{B^ 
D+TT-) = 2.9x10-3 and B{B° D*+tt-) = 2.6x10-3 0, and the height and 
the profile of dV /dq^ at large rj, or equivalently, its values at = ^ 



12 . In 



and at q^ = Mfj^ |]TT|, Md, = 1.97 GeV being the Dg meson mass 

principle, we can determine the D^*'> meson wave functions and the CKM 
matrix element \Vcb\ completely from the data fitting. 
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If there are more data from other decay modes, such as -B ^ D^*^Ds, 
fo can be fixed phenomenologically, and needs not to be specified at the 
beginning. However, these data still suffer large errors |TB|, and it is not 



practical to perform the fitting based on them. On the other hand, the data 
of S ^ p decays, B{B^ -> p~) = 8.1 x lO^^ and B{B° D*'^P~) = 
7.4 X 10~^ 0, do not give more information than B —>■ ir decays, and are not 
employed here. This is obvious from the equality of the ratios 

B{B^ ^ D^p-) ^ B{B^ ^ D+TT-) ^ 

B{BO D*+p-) ^ B{B^ D*+7i-) ^ ' ' ^ ^ 



p 



Hence, the data oi B ^ p decays just lead to the p meson decay constant 
/tt X y^8. 1/2.9 = 0.155 GeV, consistent with the currently accepted value. 
Equation (^) is a direct consequence of the vertex factorization hypothesis 
for the negligible p meson mass. Note that the Bauer-Stech-Wirbel (BSW) 
method |]TB| gives different predictions from eq. (0), which are |^ 

BjB'^D+p-) B{B' ^ D+TT-) 

Our plan is summarized as follows: 

1. Assume a set of initial values of Cd and C^*, say, Co = —2.9 GeV^ [|| 
and Cd* = —3.4 GeV^. Determine the value of Jd* from the ratio B{B^ — *■ 
D+7r^)/B{B^ D*+n-) = 1.1, which is independent of \Vcb\. 

2. Extract \Vcb\ from the magnitude of the spectrum dV /dq^ at = 
using fi)* from Step 1, and the initial Co and Co*- 

3. Determine a new value of Co from B{B^ D~^tt~) = 2.9 x 10~^, and 
a new Cd* from dV/dq^ at q^ = M^^ using fo* and \Vcb\ obtained from the 
above two steps. 

4. Go to Step 1, starting with the new initial values of Cd and Cd*- 

At last, the four parameters approach their limits after a few iterations. The 
results are 

Cd = -2.6179 GeV^ , Nd = 13.8 GeV=^ , Jd = 0.14 GeV , 
Cd* = -3.0421 GeV^ , Nd* = 14.6 GeV^ , Jd* = 0.129 GeV , 
\Vcb\ = 0.043 . (26) 

The dependence of the B and D^*^ meson wave functions on the momen- 
tum fraction x is shown in fig. 2. (pB peaks at x ~ 0.05, and (PdM peaks at 
X ~ 0.2, indicating that the heavier B meson is strongly dominated by soft 
dynamics. The profiles of 0d and (pD* are very similar, but the maximum of 
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0D* locates at a slightly smaller x compared to (pD because of the relation 
\CD*\lMf)^ > ICdI/M"^. All these features are consistent with the expecta- 
tion from the ordering of the masses, Mb ^ M^* > M^,. That the height of 
0/5 is larger than that of (po* is due to fn > fo*- Note that our prediction 
Id*/ Id = 0.92 is contrary to fo*/ Id = 1-28 appearing in the literature [|l^ . 

The parameters in eq. (|26D lead to the branching ratios B{B^ — > D^7r~) = 
2.89 X 10-3 and B{B° D*+7r^) = 2.61 x 10-^ and the spectrum dT/dq'^ 
for the semileptonic decay 5*^ D*'^i~iJ as in fig. 3, where the CLEO data 



from Ref. JlTI are also shown. It is observed that our predictions are in good 
agreement with the data at low g^, and begin to deviate above = 4 GeV^, 
the slow recoil region in which PQCD is not reliable. In order to justify the 
PQCD analysis for < 4 GeV^, or 77 > 1.3 approximately, we exhibit the 
dependence of ^+ and on the cutoff be = bi = 62 for r] = 1.3 in fig. 4. 
More than 50% of the contributions to the form factors come from the region 
with b < 0.6/A, le., ^^(l/y/Tr < 0.5. 

Variation of the six transition form factors with the velocity transfer rj is 
displayed in fig. 5. From fig. 5a, we find that the magnitudes of ^y, and 
are almost equal, with the relation = > C,Ai, and their behaviors 
are close to that of ^4., corresponding to the similarity between the profiles of 
and (pD- This similarity is also reflected by the fact that the ratio of 



to ^+ is roughly the same as fn*/ fo- Contrary to fig. 5a, ^_ and shown 
in fig. 5b increase with rj. ^_ possesses a smaller slope, and is expected to 
become negative at low rj. These features are consistent with the predictions 
from HQS combined with (9(as) and 0{1/M) corrections in ||^. 

The CKM matrix element |V"cfe| = 0.043 extracted here is a bit larger 
than recent estimations in the literature, which range from 0.035 to 0.040 
[0, |Tl], ITSf. Refer to Ref. [|l5] in which ^y, and were modeled by the 



single form factor, as indicated in eq. 



T{^) = T{l)[l-a\ij-l)+b{Ti~l) 



(27) 



with the parameters = 0.84, 6 = for a linear fit and = 0.92, b = 0.15 
for a quadratic fit to experimental data. The normalization ^(1) = rjA^i^) + 
0((A/M)^), where rjA is a perturbatively calculable quantity, takes the value 
0.93 PI, ^(1) 



^(1) 



0.89 iig 



211. In Ref. 



the 



_ _ or JP'(l) = 0.96 _ _ 

single form factor was expressed as the overlap integral of the B and D*-*^ 
meson wave functions derived from the Bethe-Salpeter equation, and was 



parametrized by a similar formula to eq. (27), 



Va 



(28) 
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for the constants rjA = 0.9942, = 1.279 and c = 0.91. The method of 
overlap integrals leads to the expression 



^iv) = exp 



2 



-2(P^ - 1 



r/ + 1 



(29) 



for the constant p2 = 1.19. In Ref. the best fit to the CLEO data [|TT 
gives the form factor 

jr(^) = 1 -0.81(77- 1) . (30) 

Comparing to our results, we find that eqs. (|27D -(^0D are in fact close to 
or larger than ^+ as shown in fig. 6. Since our form factors involved in 
the decay B — > D*ii> are smaller, the prediction of \Vcb\ is of course larger. 
If substituting ^+ for ^y, (,Ai and in eq. (^B]), we shall have \Vcb\ ~ 
0.043 X 0.9 = 0.039 extracted from the data, which then locates in the above 
range. 

We explain why our form factors are smaller than those in the literature. 
The reason is attributed to the choice of the decay constants fs = 0.12 GeV 
and /d = 0.14 GeV at the beginning of the analysis. If fs and fo increase to 
0.13 and 0.15 GeV, respectively, fc* will become 0.138 GeV because of the 
ratio fn* / Id = 0.92. Here we suppose that the shape parameters Cg, Co and 
Cn* change only slightly. Then \Vcb\ decreases to 0.037 in order to maintain 
the height of the spectrum. We have confirmed this argument by following 
Steps 1 to 4 explicitly as stated above. Therefore, the best conclusion for 
our study is that we have disentangled the task of determining \Vcb\ to the 
extent that \Vcb\ is given, in terms of and fr,, by 

/0.12 GeV\ /0.14 GeV\ , , 

for fs and fo varying around 0.12 and 0.14 GeV, respectively. Once the 
precise measurement of the decay constants fs and /d is performed, the 
CKM matrix element \Vcb\ can be fixed uniquely. 

At last, we compute the branching ratio B{B^ D^D^) inserting the 
parameters in eq. (|26|) , and compare results with the data 8.0 x 10"^ pf, if 
ignoring the errors. The expression for the decay rate is written as 

1 



X 



H 2 



(32) 
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for the CKM matrix element \Vcs\ = 1 [T^, r' = Mr) J Mb and ?7max = 
(1 + — r'^)/(2r). We obtain the Ds meson decay constant Jd^ = 0.137 
GeV, or in terms of the ratio, foslfn = 0.98. Because of the large errors 
associated with the data, we do not compare this ratio with those from the 
lattice calculation [|l^ and from QCD sum rules p2|, which are about 1.1 and 
1.2, respectively. However, our prediction is in agreement with the simple 
formula 

^ = — ^ — - ^ 0.98 33) 

for the current quark masses = 10 MeV, rUs = 150 MeV and = 1.5 
GeV. Equation ( ^3l) was obtained using general arguments from the Wigner- 
Eckart theorem and assuming that chiral symmetry is only broken by quark 
mass terms. 



4 Nonfactorizable Contributions 



In the study oi B ^ D^*^ decays, we have considered only factorizable con- 
tributions from external ly-emission diagrams as in fig. la. For the nonlep- 
tonic decay D'^n', there are also nonfactorizable contributions from 

W^-exchange diagrams as shown in fig. 7a. To justify our study, we should 
have a convincing argument that such nonfactorizable contributions are in- 
deed unimportant. A simple investigation shows that the PQCD formalism 
can be applied to fig. 7a equally well, with the following modifications: 

1. The color flow in fig. 7a differs from that in fig. la. This distinction leads 
to a factor 1/3 for nonfactorizable contributions. 

2. From the viewpoint of factorization theorems, fig. 7a is a crossing in the 
s and t channels of fig. la, excluding the color flow. That is, these two 
diagrams are similar to each other, except the interchange of the B meson 
and pion kinematic variables. Hence, the hard gluon propagator in fig. 7a is 
proportional to I/X3X2, X3 being the momentum fraction of the pion, which 
comes from the replacement of Xi by x^ in the gluon propagator l/xiX2 
associated with fig. la. Since xi and 0:3 are of order 0.05 and 0.5, respectively, 
from the B meson and pion wave functions, the interchange leads to a factor 
1/10. 

Certainly, the presence of in the hard scatterings moderates the difference. 
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Therefore, we estimate that nonfactorizable contributions from l^-exchange 
diagrams are roughly 5% of factorizable ones, and the neglect of fig. 7a is 
reasonable. 

With the parameters determined in Sect. 3, we can also study charged 
B meson decays, such as B~ D^tt~ , based on the PQCD formalism em- 
ployed here. In this case there are additional contributions from internal 
PF-emission diagrams as shown in fig. 7b. Following the similar reasoning, 
the hard scattering associated with fig. 7b, which is proportional to l/xiXs, 
is obtained by interchanging the kinematical variables of the D meson and 
of the pion. This interchange gives a factor 0.2/0.5 = 2/5. We estimate that 
contributions from internal Pi^-emission diagrams, combined with the color- 
suppressing factor 1/3, are roughly 15% of factorizable ones, which are of 
course sizable. This conclusion is consistent with predictions from the BSW 
method [^]. We shall discuss these subjects in details in a seperate work [P^]. 



5 Conclusion 



In this paper we have fixed the D^*^ meson wave function using the exper- 
imental inputs from the nonleptonic decays B — > D^*\, and evaluated the 
spectrum of the semileptonic decay B D*lv at large velocity transfer. 



IK 



cb\ 



0.043 is extracted. The form 



from which the CKM matrix element 
factors involved in i? — > D^*'^ transitions are obtained, and compared to the 
predictions from HQS combined with 0{as) and 0{1/M) corrections 
from the overlap integrals of heavy meson wave functions P, |10[, and from 



the data fitting [|Tl|]. The value of \Vch\ extracted here is larger than those 
in the literature [IC, 11, 18|, and the reason is that we have adopted the 
smaller decay constants and fr). A precise measurement of fs and fo in 
the future will remove this ambiguity. 

We emphasize that the behaviors of all the transition form factors are 
derived from the single D*^*^ meson wave function that is fixed at specific 
kinematic points, without resorting to a model for each form factor, such 
as the algebraic forms employed in the pole forms in |]T6|, ^ and the 
exponential forms in p6|. Note that all the above model form factors are 
larger than ^_|_ presented in this work. Hence, they lead to smaller values of 
\Vcb\-i ranging from 0.032 to 0.038. Our formalism, based on the parameters 
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in eq. p^), can be used to study B meson decays, especially the nonleptonic 
cases with nonfactorizable contributions, in a less ambiguous way . 

We thank H.Y. Cheng and P. Kroll for useful discussions. This work 
was supported by National Science council of R.O.C. under the Grant No. 
NSC-85-2112-M-194-009. 



Appendix 



In this appendix we present the explicit expression of the exponent s{k, h) 
appearing in eq. (|^). The full expression, instead of the first six terms [^], 
is adopted in this paper. It is given, in terms of the variables. 



by 



q = \n{k/K), 6 = ln(l/6A) 



AW 



2Pi 



gin ( S 



' 8/5? 

" 16/3^ 
v4(i)/32 q - b 



In 



,27-r 



In 



4/?2 

:i) 

ln(2g) + 1 ln(26) + 1 
q b 



1-1 



In 



ln(2g) + 1 ln(26) + 1 



21n(2g) + 3 2 ln(26) + 3 



1728/3f 

A^q_-_ 
432/5f 63 



2 ln(26) + 1 



(34) 



62 

181n2(2g) + 301n(2g) + 19 18 W(26) + 30 ln(26) + 19" 



q^ 



52 



91n2(26) + 61n(26) + 2 



(35) 
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The above coefficients (3i and A*^*-* are 




/32 



24 

67 TT^ 10 
Y ~ y ~ 27 



153 - 19ra 



8 /e'^-^\ 



where 7^; is the Euler constant. 

Note that s is defined for q > b, and set to zero for q < b. As a similar 



treatment, the complete Sudakov factor exp(— 5*) is set to unity, if exp(— S*) > 
1, in the numerical analysis. This corresponds to a truncation at large kr, 
which spoils the on-shell requirement for the light valence quarks. The quark 
lines with large /c^ should be absorbed into the hard scattering amplitude, 
instead of the wave functions. An explicit examination shows that the partial 
expression including only the first six terms gives predictions smaller than 
those from the full expression by less than 5% [@]. 
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Figure Captions 

Fig. 1. (a) Factorization oi B ^ D^*^ transitions, (b) 0{as) corrections to 
wave functions. 

Fig. 2. Dependence of the B and D^*^ meson wave functions on the mo- 
mentum fraction x. 

Fig. 3. The spectrum dV /dq^ of the semileptonic decay B^ D*^i~i?. 

Fig. 4. Dependence of ^+ and C,Ai on the cutoff be for r] = 1.3. 

Fig. 5. Dependence of (a) ^v, and and of (b) ^_ and on r]. 

Fig. 6. Dependence on rj of (1) eq. (|9D, (2) eq. (HD, (3) eq. (ID, (4) 
eq. ( P7D for the linear fit with J-'(l) = 0.93, and (5) eq. (^) for the quadratic 
fit with = 0.93. The curves of ^+ and C,v are also shown. 

Fig. 7. (a) ly-exchange and (b) internal ly-emission 0{as) diagrams. 
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